THE HEAVY PARTICLE HAZARD - WHAT PHYSICAL DATA ARE NEEDED ? 
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It is impossible at present to evaluate adequately the radiation hazard 
from heavy galactic cosmic rays to astronauts on extended missions. Both 
physical and biological data are sorely needed for the proper assessment 
of the risk on missions lasting a year or longer outside the earth's magnetosphere . 
The physical data required fall into three main categories: 

1. Spectral Characteristics. From the relative abundance of the highZ nuclides 
and the limited data already obtained at solar minimum in the low energy 
region, it is clear that the flux of very high LET particles (above 2 x 1 0 3 

MeV cm / g) will be predominantly low energy ions « 500 MeV/nucleon) 
in the iron group (26 < Z < 28) . Thus, it is important to determine the shape 
of the low energy portion of the spectrum during the next solar minimum 
(~1976), when galactic fluxes reach their maximum values. 

2, Nuclear Interaction Parameters. Although ionization will be the main 
energy loss process for the heavy ions, a considerable fraction will also 
undergo nuclear collisions, fragmenting into particles of lower charge but 
perhaps greater range than the primary particle. Because the mean free 
path for interaction (-vl5 g/cm ) is on the order of the probable shielding 
available for critical body organs (the eye, brain), it is necessary to know 

the secondary production characteristics in order to make a reasonable estimate 
of particle fluxes within the body organs in question. Results of a computer 
code are presented in which two different assumptions for the fragmentation 

parameters and their Z dependence (experimentally unknown at present) have 
been made. 
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3. Track Structure. The spatial distribution of ionization surrounding a high 

energy, high 2 ion plays an important role in determing the ultimate 

biological damage. Knowledge of track structure will help to predict 

biological effects and may be crucial in the all-important extrapolation 

from animals to man. A rough theory for calculating the dependence of the 

energy density in water as a function of distance from the track trajectory 

will be reviewed. Limited biological data at much lower energies suggest 

2 2 

that the height of the shoulder of this curve (which is proportional to 2 / B , 
not total LET) may be closely correlated to biological damage. 

Data in the first category above must be obtained by spacecraft flying 
outside the magnetosphere during solar minimum and instrumented to measure 
the low energy portion of the high Z galactic component. The data in the 
second and third categories can be most reliably and readily obtained in 
earth-based laboratories (along with the necessary biological data). The 
requirement here is an accelerator of ions in the iron group up to the region 
of 300-500 MeV/ nucleon. 


A convincing case has been made that on re- 
cent Apollo missions, astronauts have experienced 
within their eyes, and probably via direct retinal 
excitation, the effects of the passage of high ener- 
gy heavy galactic cosmic rays (ref. 1) . It remains 
to be determined the extent of any lasting, dele- 
terious effects produced by such ions over a con- 
tinuous extended exposure. At the present time, 
there are no data on the biological effects of such 
high LET penetrating radiation on man or animals. 
The reason for this, of course, is that no acceler- 
ator exists at present which can produce beams of 
heavy ions with sufficient penetration to attempt 
the type of biological experiments necessary for 
the proper evaluation of this hazard. There is 


much to learn before we ca n be sure that astronauts 
on very long missions, that is, missions of a year 
or longer, will be safe from the continuous bom- 
bardment of these high LET particles. 

In this paper we will not deal with the biolo- 
gical research that must be done to study the long 
term effects of such radiation, but insteadwill dis 
cuss the physical data necessary for an adequate 
evaluation of the problem. 

The research falls into three main categories, 
each of which is interesting in its own right. These 
are presented in Table I. 

First, there is the radiation environment itself 
what is the flux of these high energy, heavily 
ionizing ions; how many might be expected to hit 
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PHYSICAL DATA REQUIRED 


RADIATION ENVIRONMENT OF HIGH ENERGY HEAVY IONS 

Fluxes 

Energy spectra (emphasis on energies ( 500 MeV/nucleon) 

Solar modulation 

NUCLEAR INTERACTION PROPERTIES 

Cross sections in aluminium and water (tissue) 

Fragmentation parameters 

TRACK STRUCTURE 

Radial energy density 

Track core size 

Correlation of structure to biological effects 

a spacecraft on, say, a Mars mission? Of particular 
interest is the energyregion just below 500MeV/ 
nucleon, where the penetration is still sufficient 
to be of concern and where the ionization (that is 
the LET) will also be high.* Solar modulation is 
important since these lower energy particles are af- 
fected more than those of higher energy by the mod- 
ulation mechanism. Long manned missions may be 
planned for solar inactive years to decrease the 
probability of encountering giant solar-particle 
events. If so, the galactic fluxes encountered wi 1 1 
be near maximum . 

Secondly, we must study the nuclear interaction 
properties of the particles involved in the materials 
of tlie spacecraft and the bodies of the astronauts 
themselves. No data are presently available on the 
cross sections of heavy ions in the 1 GeV/ nucleon 
range in aluminum and water or tissue. The total 
non-elastic cross sections can be estimated with 
fair accuracy, but we mustalsoknow the fragmen- 
tation parameters associated with the various par- 
ticle reactions. The fragmentation parameter for a 
given secondary particle type is either the proba- 
bility of emission of the secondary, if less than unity, 
or the average number of secondaries per interac- 
tion, if greater than unity. Also, how is the energy 
shared as the two nuclei break up ? 

* For orientation, a 500 MeV/nucleon iron nucleus ha 
of 2000 MeV cm 2 /g (200 keV/jj) . 


Does the resulting fragment continue with the same 
velocity, oristhere significant change in the energy 
per nucleon of the emerging heavy particle ? 

These questions must be answered before a clear 
idea can emerge of the penetration characteristics 
of the heavy ions. Only then will it be possible to 
make accurate predictions of the fluxes of high Z 
ions at organs which may be critical in that they 
might accumulate enough damage to impair normal 
function . 

Finally, there is the field of track structure. 

Here the interest is in the spatial distribution of 
energy deposition by the heavy ions as they slow 
down in tissue. What is the extent of the core of 
such an ion ? What is its “radius of destruction" ? 

Is the important physical parameter its LET or dE/dx, 
or does the diffusive character of the "electron 
penumbra" decrease its biological effectiveness to 
less than that of a lower energy ion at the same 
dE/ dx? These are all questions which bear directly 
on the biological effects of these ions. Since labora- 
tory experiments may be performed with animals 
which have cells or cell nuclei with different sizes 
than the equivalent cells in man, it may be impor- 
tant to know as much as possible about the physical 
distribution of energy density around a heavy ion 
track in order to extrapolate results from animals 
to man. 

We now review briefly some of the things we 
already know in these areas. 

Integral Number-LET Spectra 

Considerable data have been obtained from satel- 
lites and balloons on the fluxes and energy spectra 
of heavy galactic particles during the last solar 
minimum period; that is, during the period of maxi- 
mum galactic fluxes - 1964-66 (ref. 2). 


2 
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In figure 1 we show the integral spectra not as 

a function of the energy or energy per nucleon as 

is usually presented/ but as a function of the dE/dx 

or LET of the particles/ which has more biological 

relevance. We note that the particles in the iron 

group (26<Z£28) dominate the spectra at very 

2 

high LET (>1000 MeV cm /g) and also that the 

data is scanty - the solid curve stops at 5000 MeV 

2 

cm /g corresponding to an energy of 100 MeV/ 
nucleon. The dashed lines at high LET for each 
component are just guesses which are consistent 
with the low energy spectral behavior of the lighter 
components we know more about - protons and 
alpha particles. The dotted line at high LET was 



Figure 1 • Integral number-LET spectrum for the gal- 
actic cosmic rays behind no shielding at solar mini- 
mum. The contribution of the various Z groups are 
also shown. The dotted line was obtained assuming a 
sharp low-energy cutoff of 100 MeV/nucleon to the 
iron spectrum. The dashed lines represent regions where 
the data are lacking. 


calculated assuming a sharp low-energy cut off in 
the iron spectrum at 100 MeV/nucleon. We con- 
clude that more data are needed on the spectra of 
the very heavy, i.e. iron group ions, in the energy 

region just below 300 MeV/nucleon (residua I range 
2 

in water: 4-5 g/cm ) during the next solar cycle 
minimum (1975-77). Since these particles are 
strongly affected by solar modulation mechanisms, 
it would also be of interest to determine when and 
to what extent the low energy fluxes of these heavy 
ions decrease during the cycle. 

Nuclear Fragmentation 

For high energy heavy ion-nuc leus interactions, 
we have only cosmic ray data to give us any feel- 
tng f or the cross sections and types of secondaries 
involved. In figure 2 we show a compilation of 

number of events in 
parentheses 



Figure 2. Compilation of experimental data for frag- 
mentation parameters of heavy ions on the light ele- 
ments (excluding hydrogen) in nuclear emulsibn. The 
data come from rather old cosmic ray studies on bal- 
loons by the Bristol group (refs. 3 and 4). 
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data gathered from nuclear photographic emulsion 
on the fragmentation parameters in the light ele- 
ments of the gelatine/ i.e. carbon, nitrogen and 
oxygen. This was all the emulsion data in the lite- 
rature a few years ago and was taken from the old 
compilations of the Bristol group (refs. 3 and 4). 
Now recent balloon flights could undoubtedly yield 
more data of relevance here, but the heavy ele- 
ments in the emulsion (silver and bromine) will al- 
ways give a certain amount of contamination trouble 
and considerable scanning would be required to get 
good statistics on each incident ion. Of lesser im- 
portance is the fragmentation in materials of which 
the spacecraft itself is made, for instance aluminum. 

In the Space Physics Group at Boeing, a compu- 
ter code has been developed to calculate the result- 
ant fluxes and dose behind various shielding thick- 
nesses for any incident galactic particle spectrum 
or combination of spectra. The results and details 
of the calculation have been presented earlier in 
this symposium (ref. 5) and will not be repeated here, 
except to show the two sets of fragmentation para- 
meters selected and the results of the dose rate cal- 
culations using the different sets. Figure 3 shows 
the two sets used. In set I, a constant probability 
for emission as a function of secondary Z is assumed, 
such that total charge is conserved. In set II, we 
try to reflect the possibility that fragments with Z 
close to two and close to the Z of the primary par- 
ticle might be more probable than intermediate Z 
secondaries. In figure 4, the results of the dose rate 
as a function of depth for the iron group spectrum 
using the two sets are shown as the upper solid and 
dashed lines. There is a factor of two difference at 
2 

60 g/cm of water. Also shown are representative 
curves of the secondary dose rates from the 
product secondary fragments. These are just repre- 
sentative curves since there are actually twenty- 
five such secondary curves plus the uncollided iron 
curve contributing to the total. 


Set 1 



secondary Z, j 


Set 2 



Figure 3. Two sets of fragmentation parameter second- 
ary Z dependence assumed in the Boeing computer pro- 
gram mentioned in the text. Set I assumes a constant 
dependence on secondary 2 plus charge conservation. 
Set II reflects the possibility that secondary ions with 
Z just less than the primary Z and secondary ions with 
Z close to that of a particles might be more probable 
than intermediate Z secondaries. 
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^Y r ) as a function of distance from the track tra- 
jectory. The integral of 2vr ft. (r) over distance 
from the trajectory yields the average dE/dx or 
LET : 


LET = 2ir 


poo 

j{V r ) rdr 


The distance r is the perpendicular distance from 
the track trajectory. It is easier to study the quan- 
t'ty ( r ) large distances, say greater than 50- 
100&, than at small distances. At large distances, 
the only electrons present have come from the so- 
called "close collisions" and can be treated with 
the familiar Rutherford expression for free electrons. 
Two rather crude theories (refs. 6 and 7), which 

will not be discussed here/ yield a dependence of 

P 2 2 

j (r) going as Z /ft of the particle in question 

Figure 4. Dose rates as a function of depth from galactic and faMing Q5 a functfon of r as ! / r 2 up to d!stances 
cosmic rays of the iron group including secondary pro- wh!ch are sMM sma || com p are d with the range of 
duction as calculated by the Boeing code. Total dose electrons with maximum transferable energy. As we 
rates as well as representative values from the second- approach this range, of course, the function falls 
ary contributions are presented. The solid linds result off much more rapid | y . Q ne estimation of the depen- 
from the Set I fragmentation parameters and the dashed dence of ^ on r for an iron nuc | eus of energy 
result from the Set II fragmentation parameters presentedsoo MeV/nucleon is shown in figure 5. Here the 
in Figure 3. Note a difference of a factor of two in the ordinate is the "dose" in units of 100 ergs/g. 


total dose rates at 60 g/crn- of water. 

We conclude that a better knowledge of the 
fragmentation parameters for such heavy Ion break- 
ups Is essential. At the same time/ of course/ this 
will provide us with more insight into nuclear 
structure and the way energy is shared in high ener- 
gy nucleus-nucleus collisions/ an area of interest 
to the nuclear physicist. 

Track Structure 

The physical distribution of energy deposition 
(ionization density) along the trajectory of a char- 
ged particle has been an area of increasing interest 
in recent years. One way of approaching this pro- 
blem is by defining an average energy density 


Dose is in quotes here because it is not a dose in 
the strict sense of the word. It can be considered 
an average dose over distances along the track large 
enough so that dE/dx has meaning/ that is/ long 
compared to distances where the fluctuation of 
energy loss plays an important role. The dependence 
of this quantity at small distances is an open question 
at present. It is assumed to rise steeply because we 
know that the total integral must equal the total 
dE/dx/ and subtraction of the tail contribution from 
the total leaves a large portion of the energy to 
be deposited in the core, at least one-half the ener- 
gy. Hence the core is undoubtedly very important 
in producing biological effects. This nonuniformity 
of energy deposition implies that the average dose 
deposited by such particles or their total LET may 
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Figure 5. Estimation of the dependence of the energy 
density on distance from the particle trajectory of a 
500 MeV/nucleon iron ion. The ordinate is not a true 
dose but can be considered an average dose over dis- 
tances along the track long compared to distances 
where fluctuations of energy loss are important* 


have little relevance to the evaluation of the hazard 

from these high-energy heavy ions. 

There is evidence at much lower energy that the 

2 2 

quantity Z /ft is a more relevant quantity than 
LEToo to ^ escr '^ e biological effects (ref. 8). It is 
well known that hypoxic cells are more radioresis- 
tant than oxygenated cells. The Oxygen Enhance- 


ment Ratio (OER) is a measure of this difference 
in radiosensitivity for a particular type of radiation. 
An OER of 2.7 as measured in human kidney cells 
for X-ray irradiation/ for instance/ means that 2.7 
times the dose is needed to produce the same effect 
on hypoxic cells as on well-oxygenated cells. If 
for a different type of radiation/ the difference in 
radiosensitivity is less/ the OER is lower. For an 
OER of unity/ there is no difference in the radio- 
sensitivity of oxygenated and hypoxic cells. 

In general/ the OER is closer to unity for radiations 
of high LET than for those of low LET. This is shown 
for two sets of experimental data in fig. 6a for the 
system of human kidney cells in tissue culture. The 
two sets of data differ in that the solid squareswere 
obtained with high-energy heavy ions (•* 7 MeV/ 
nucleon) all at the same velocity (ref.9),and the 
open circles were obtained with alpha-particles 
at different velocities (ref. 10). We note that there 

is no overlap of the data. If 7 however, we plot 

2 2 

the data against Z* / ft instead of against LET^, 
we obtain the results in figure 6b. Here Z* is the 
charge of the particle, taking into account charge 
pick-up of the ion as it slows down. The data over- 
lap to a much greater extent and are consistent 
with lying on the same curve. Thus we conclude 

that at least for this biological system and end- 

2 2 

point, Z* / ft is a more relevant physical quantity 
to plot the biological results against. This corre- 
lation should be checked with other endpoints and, 
in particular, should be checked with high-energy 
heavy ions of importance in the galactic cosmic 
rays, where the "electron penumbra" is greater 
in extent, thus diffusing part of the energy deposited 

over a greater distance. The "close in" shoulder 

2 2 

of the curve, as seen in fig. 5, rises as Z / ft , 

2 

so for high Z the shoulder will also be large. 

Core effects are extremely important because of 
the very high energy density in the core and it 
is even possible that new inactivation processes 
will come into play which are entirely absent at 
lower energy densities. 
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Figure 6. a. OER of human kidney ceils ( in vitro ) from 

heavy ions at ~7 MeV/nucleon (solid squares) and from 

low energy alpha particles (open circles) plotted against 

LET • b. Same data as in a. plotted instead against 
oo 2 

Z*2/B , where 2* is the effective charge of the ion 
(including charge pick-up) and ft is the ratio of the 
velocity of the particle to that of light. 


Conclusion 

We have mentioned several areas where physical 
data are needed in order to evaluate adequately 
the potential hazard of high-energy heavy ions in 
the galactic cosmic radiation to long manned 
missions outside our magnetosphere. The first/ more 
data on the lower energy portion of the differential 
energy spectrum of the iron group/ can be obtained 
in appropriately instrumented satellites during the 
next period of solar inactivity 1975-1977. The 
other two areas, nuclear interaction parameters in 
tissue and aluminum and the structure of the energy 
deposition along the trajectory of the tracks, appear 
to require an accelerator which accelerates ions of 
the iron group up to at least 300-500 MeV/nucleon . 
This information, along with the necessary biological 
information concerning the extent of damage and 
functional impairment to various critical organs, is 
necessary before we can confidently send manned 
expeditions to explore the near planets or establish 
bases for extended exploration of the moon. 
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